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Abstract. We present a comprehensive study of the Magellanic Cloud planetary nebula SMP61 and of its nucleus, a Wolf- 
Rayet type star classified [WC 5-6]. The observational material consists of HST STIS spectroscopy and imaging, together with 
optical and UV spectroscopic data collected from the literature and infrared fluxes measured by IRAS. We have performed a 
detailed spectral analysis of the central star, using the Potsdam code for expanding atmospheres in non-LTE. For the central star 
we determine the following parameters: = IO^'^Lq, R^, = OAIRq, = 87.5 kK, M = 10"* Mq yr"', v„ = 1400 kms"', 
and a clumping factor of D = 4. The elemental abundances by mass are Xhc = 0.45, Xq = 0.52, < 5 10"', Xq = 0.03, 
and Xpc < 1 lO""*. The fluxes from the model stellar atmosphere were used to compute photoionization models of the nebula. 
All the available observations, within their error bars, were used to constrain these models. We find that the ionizing fluxes 
predicted by the stellar model are basically consistent with the fluxes needed by the photoionization model to reproduce the 
nebular emission, within the error margins. However, there are indications that the stellar model overestimates the number and 
hardness of Lyman continuum photons. The photoionization models imply a clumped density structure of the nebular material. 
The observed C m] /11909/C ii A4267 line ratio implies the existence of carbon-rich clumps in the nebula. Such clumps are likely 
produced by stellar wind ejecta, possibly mixed with the nebular material. We discuss our results with regard to the stellar and 
nebular post-AGB evolution. The observed Fe-deficiency for the central star indicates that the material which is now visible 
on the stellar surface has been exposed to s-process nucleosynthesis during previous thermal pulses. The absence of nitrogen 
allows to set an upper limit to the remaining H-envelope mass after a possible AGB final thermal pulse. Finally, we infer from 
the total amount of carbon detected in the nebula that the strong [WC] mass-loss may have been active only for a limited period 
during the post-AGB evolution. 

Key words, stars: Wolf-Rayet - stars: atmospheres - stars: mass-loss - ISM: abundances - ISM: planetary nebulae: individual: 
SMP61 - planetary nebulae (general): ISM: abundances 



1. Introduction 

Only a few studies have been devoted so far to a consistent 
modelling of a planetary nebula and of its central star Such 
studies are useful to get a better insight into the relation be- 
tween the nebula and its progenitor Another, very important 
aspect is that this is the only way to test model atmosphere 
predictions in the Lyman continuum and thus to validate the 
model atmospheres. The works of Rauch, Koppen & Werner 
(1994, 1996), Pefia et al. (1998), De JVIarco & Crowther (1998, 
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1999), De Marco et al. (2001) are examples of such studies, 
while Crowther et al. (1999) have performed a similar study 
on a Population I Wolf-Rayet ring nebula. Such investigations 
are particularly important in the case of planetary nebulae with 
Wolf-Rayet type central stars (which represent about 10% of 
all planetary nebulae), since recent work (e.g. Gorny & Tylenda 
2000, De Marco & Soker 2002) have completely changed pre- 
vious views on the evolutionary status of these objects. 

In general, Wolf-Rayet central stars of PNe belong to the 
[WC] sequence. In our galaxy most of these objects have been 
classified as [WC-early] or [WC-late] types, with only few ob- 
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jects of intermediate types (Tylenda, Acker & Stenholm 1993). 
In the Magellanic Clouds, the WR central stars are also of 
[WC] type (except for the extraordinary central star of LMC- 
N66, e.g. Pena et al. 1997b), but in this case they have been 
classified in the intermediate [WC] types. Pena et al. (1997a) 
suggested that such a difference might be a consequence of 
the differences in metallicity between the Milky Way and the 
Magellanic Clouds. In any case, [WC] central stars show spec- 
tral features identical to those of massive WC stars but at much 
lower luminosity, and they can be analyzed with the tools de- 
veloped for massive WR stars (e.g. Hamann 1997). 

In the present paper, we concentrate on the planetary nebula 
SMP61 (also known as N203, WS 24 and LMl-37) which is 
among the brightest planetary nebulae in the Large Magellanic 
Cloud. This is a good case for a detailed st udy: The o bject is at 
a known distance modulus of 18.50 mag ( Benedict et al..2002) . 
equivalent to 50. 1 kpc. The central star is one of the brightest 
[WC] central stars in the LMC and therefore relatively easy to 
observe. It is of [WC5-6] type (Monk, Barlow & Clegg 1988; 
Pena, Ruiz & Torres-Peimbert 1997a), therefore the nebular 
spectrum contains many lines of various excitation levels that 
allow refined diagnostics. In addition, the nebula appears to be 
spherical and integrated spectra are already available (Pena et 
al. 1997a). 

Using HST STIS, we have secured high signal-to-noise 
spectra of SMP61 in a wide spectral range. This provided 
strong constraints for our modelling of the central star atmo- 
sphere. The best fit model atmosphere was then used as an input 
to build a photoionization model for the nebula. The observa- 
tions are described in Sect. 2, the stellar modelling technique is 
presented in Sect. 3, the spectral analysis of the central star of 
SMP61 is discussed in Sect. 4. The nebular modelling strategy 
is exposed in Sect. 5 and the nebular model fitting of SMP61 
is presented in Sect. 6. The implications of our modelling are 
discussed in Sect. 7, and the main points of this study are sum- 
marized in Sect. 8. 

2. New Observations and Previous Data 

2.1. Observations 

HST STIS spectroscopic data were obtained on 1998/12/30 
(HST Cycle 7, program ID 7303). MAMA and CCD detec- 
tors were employed to cover the broadest spectral range possi- 
ble. Prior to the spectroscopic observations, STIS images were 
obtained to acquire the target. The log of the observations is 
presented in Table 1 . 

Direct images with HST WFPC were searched in the 
archives. Apart from our own STIS images for target acqui- 
sition (see Table 1), we retrieved the FOC image W1ID0601T, 
obtained on 1993/11/08, with exposure time of 3 min through 
the filter F502N which isolated the strong [O iii] A5007 line 
(this image is shown in Vassiliadis et al. 1998b). We decided 
to use our target acquisition images to analyze the nebular mor- 
phology and surface brightness because they are better exposed 
and because the Vassiliadis et al. images were obtained before 
the COSTAR mission. Our images are obtained with a filter 
including the brightest emission lines: Ha, H/3 and [O m] lines. 



We also considered ground-based spectrophotometric data 
for SMP61 obtained by Monk, Barlow, & Clegg (1988), 
Meatheringham & Dopita (1991) and Pena et al. (1997a). Other 
available UV spectroscopic data are lUE data given in Pena et 
al. (1997) and HST FOS data from 1 175 to 4800 A published 
by Vassiliadis et al. (1998a). 

All the available spectroscopic data are presented in 
Tabled which lists the reddening corrected fluxes for the ob- 
served lines (in units of H/?=100) obtained by the different au- 
thors. We also list the total H/3 flux observed in the slit, F(HjS), 
the logarithmic extinction at Hj3, C(Hj3), and the slit dimen- 
sions, when available. All the spectra mentioned above have 
good signal-to-noise, and the differences among observations 
by different authors are, in general, small. Note that observa- 
tions by Pena et al. (1997a) were made both with a slit of 2"and 
a slit of 10"and, as noted by these authors, the line ratios are 
very similar in both slits. The H/3 flux reported in column 2 of 
Table|2lis the one obtained with the 10"slit. 

SMP 61 is one of the twelve planetary nebulae in the LMC 
detected by IRAS. The fluxes retrieved from the IRAS data 
base are: F(12/im) = 0.08 Jy (flux density quality = 2) F(25/vm) 
= 0.13 Jy (flux density quality = 3), F(60yL(m)<0.16 Jy (flux 
density quality =1). 

2.2. Preliminary analysis of the observational data 

In HST images the nebula appears perfectly spherical and 
shows no small-scale structure. The surface flux distribution, 
as derived from the image o57n02iuq, is shown in Fig .[2 The 
ionizing star dominates the emission in the central pixels. The 
nebula seems to have an inner radius of about 0.065", show- 
ing the maximum surface brightness there. The surface bright- 
ness declines outwards and drops below 1/10 of its maximum 
at 0.26". 




a. 5 5 7.5 



Fig. 1. Radial distribution of the surface brightness. The STIS 
image 57n02iuq, obtained for target acquisition, was employed 
to prepare this graphic. The radius is measured in pixels. Each 
pixel is equivalent to 0.05". 

The extinction of the nebula is modest. Throughout the pa- 
per, we use the value C(H/3) = 0.18 (corresponding to E^-v - 
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Image ID 


Date 


A Range 
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^ 1 /I 9/QzL 
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i U o VV 






1 1 S9 1 Q7S 




Large ap., low res. 




UJ / IIkJjLKJ IXj 








FT IV-MAMA 


HST STIS G140L 


o57n02020 


07/01/99 


1124-1730 


22 


FUV-MAMA 


HST STIS G230L 


o57n02030 


7/01/99 


1580-3160 


25.5 


NUV-MAMA 


HST STIS G230L 


o57n02040 


7/01/99 


1580-3160 


54 


NUV-MAMA 


HST STIS G430L 


o57n03010 


30/12/98 


2900-5710 


37.86 


CCD 


HST STIS image 


o57n02iuq 


06/01/99 


5852 


3.3 


CCD, acquisition im. 


HST STIS image 


o57n03bqq 


30/12/98 


5852 


3.3 


CCD, acquisition im. 



(a) CTIO 4m with Reticon detector and grating KPGL2 

Also in HST archives FOS spectra Y2N30302T, Y2N30303T, Y2N30304T, 1 150 - 4800 A 
reported by Vassiliadis et al. (1998b) 



0.122 mag derived by Pena et al. (1997a) from the observed 
Balmer decrement and we adopt the Seaton (1979) reddening 
law. 

The STIS observations were essentially obtained to provide 
constraints on the central star. However, the 0.2"x52"STIS slit 
includes a significant fraction of the nebula so that both stel- 
lar and nebular emission is seen. The nebular emission must of 
course be subtracted from the stellar emission in order to per- 
form the stellar analysis (see section 4.1). Most of the emission 
lines detected in the UV are of stellar origin, except for the in- 
tense Cm] 1909 A, Cii 2321,26 A, and [On] 2470 A. Many 
bright nebular lines were detected in the optical region. The 
dereddened fluxes of all the nebular lines measured in the STIS 
calibrated spectra are presented in the first column of Tabled 

The STIS nebular lines together with the data presented 
by Pena et al. (1997a, column 2 in Tabled) have been used 
to derive first estimates of the physical conditions and of the 
chemical composition of the ionized gas. The atomic data are 
the same as used in the nebular modelling (see Sect. 5). The 
density derived from [S n] /16731/6717 is 5400 cm"-', the tem- 
perature derived from [O iii] /14363/5007 is 10300 K, and the 
one derived from [Nii] .15755/6584 is 11400K. These values 
were used to derive the ionic abundances, and the ionization 
correction factors from Kingsburgh & Barlow (1994) led to the 
following total abundance ratios : H : He : C : N : O : Ne : S 
: Ar = 1: .104 : 7.35 10"'* : 3.95 10"^ : 2.65 10"^ : 4.26 10"^ : 
6.0 10"* : 8.0 10"^. These abundances are employed in Sect. 6 
as starting values for the photoionization modelling. 

Regarding the IRAS data, Zijlstra et al. (1994) commented 
that the IRAS colors F(12)/F(25)~0.62 and F(25)/F(60)<0.81 
of SMP 61 are very blue and typical of a young nebula. 



3. Stellar modelling technique 

3.1. Model atmospheres 

The analysis of hot stars with strong winds like the [WC]- 
type central star of SMP 61 requires the application of 
model atmospheres that account for complex model atoms 
in non-LTE. Recent investigations revealed the necessity to 



include iron line-blanketing as well as the elfect of den- 
sity inhomogeneities (clumping) for a reliable determina- 
tion of their stella r parameters and emergent flux distri- 
butions fa amann & Koesterkd Il998t iHiUier & MiUeil 1 19991: 
ICrowther et alJ2002HGrafener et alJ2002i) . 

In the present work, we utilize the Potsdam code 
for e xpa nding atmosp heres in no n-LTE ((Koe sterke etal 



1994; 


Koesterke & HamannI 1 19951 iLeuenhasen et al 


1996; 


Hamann & Koesterke 


Il998l iKoesterke et al. 




2002; 



iGrafener et alJ l2002h . This code calculates the radiative 
transfer in the co-moving frame of reference for a spherically 
symmetric, stationary outflow. The atomic populations, the 
electron density, and the electron temperature throughout the 
extended atmosphere are determined from the equations of 
statistical- and radiative equilibrium. The system of statistical 
equations is solved in line with the radiation transport by 
application of the ALI forma lism (accelerated lambda itera- 
tion, see Hamanr] [l98lll986l) . whereas radiative equilibrium 
is obtained by a tempe rature c orrection pr ocedure whi ch is 
based on the method o f'Unsolc^ (F l955i) an d ILucvI (1 1964 (see 
Hamann & Griifener 2003; Mihalad l 19781 p.l74). The statis- 
tical equations are solved for complex model atoms of He, C, 
O, Si, and the iron group, taking advantage of the concept of 
super levels fo r the inclusion of m illions of iron-group line 
transitions (see iGrafener et alJ2002h . Density inhomogeneities 
are accounted for by the assumption of small-scale clumps 
with a constant volume filling factor fv ( Haman n & Koesterke. 
Il998l). 



3.2. Model parameters 

The model atmospheres are specified by the luminosity and 
radius of the stellar core, by the chemical composition of the 
envelope, and its density- and velocity structure. The basic pa- 
rameters are: The stellar core radius at Rosseland optical 
depth tr = 20, the stellar temperature (related to the lumi- 
nosity Li, via the Stefan Boltzmann law), the chemical com- 
position (given by mass fractions Xhc, ^c, ^o, ^si, and Xp^), 
the mass-loss rate M, the terminal velocity Voo, and the "clump- 
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Fig. 2. Subtraction of the nebular con- 
tinuum. Top panel: The observed flux 
distribution around the Balmer jump be- 
fore and after subtraction of the nebular 
continuum. Bottom panel: The reddened 
theoretical continuum as utilized for the 
nebular subtraction is plotted in black. It 
is compared to the nebular flux (grey) as 
estimated from the observation, by cor- 
recting for the stellar contribution (see 
text). Although the latter may still in- 
clude a small contribution from the cen- 
tral star, the theoretical flux distribution 
is clearly confirmed. 



ing factor" D - I /fv giving the density enhancement of the 
clumped matter with respect to the mean wind density. For the 
dependence of the velocity on radius, a relation of the form 
v(r) = Voo (l-Ro/r) is assumed. This outer velocity law is 
augmented by an exponentially decreasing density distribution 
in the hydrostatic domain. Ro is suitably determined to connect 
both domains smoothly. 

The parameters M, Vco, and D are connected by the 
"transformed radius" Rt oc 7?*(voo/ VdM)^^\ Models with 
the same Rt show a lmost identical line equivalent widths 
jSchmutz et aljfl989l) . This invariance holds because the line 
emi ssion in WR stars is dom inated by recombination processes 
(see iHamann & Koesterkell998) . 



4. Spectral analysis of the central star 

In the present section, we perform a detailed spectral analysis 
of the central star of SMP61 based on HST UV + optical ob- 
servations. Owing to the reliable spectrophotometry of HST in 
combination with the known distance and the moderate inter- 
stellar reddening towards the LMC, the observations provide 
the absolute stellar flux distribution in the wavelength range 
from 1150 to 6000A. The model comparisons in the present 
section are therefore performed with respect to the absolute 
flux, c orrected for interstellar extinction applying the standard 
law of lSeatonl(tl979.) with Eb-v =0.122 mag. 

4.1. Subtraction of the nebular continuum 

Due to the small angular diameter of the nebula (0.7"), the HST 
spectra (0.2" slit width) include a significant part of the to- 
tal nebular flux. Consequently, the observations clearly show 
nebular features superposed to the stellar spectrum. The most 
striking of these is the Balmer jump at 3646 A in emission, in- 
dicating the dominance of the nebular continuum in the UV 
to blue part of the observation. A reliable subtraction of this 



continuum emission is essential for the spectral analysis of the 
central star 

A distinction between both contributions is in principle 
possible, because the nebula is spatially extended in contrast to 
the point-like star. To estimate the contribution of the nebula, 
two spectra with apertures of 1 .2" and 0.25" are extracted from 
the UV part of the HST data. The difference between these two 
spectra is dominated by nebular emission. After scaling the dif- 
ference spectrum to the flux of the nebular C ii line at 2327 A, 
a coarse estimate for the pure nebular flux is obtained. 

In Fig.|2] this spectrum is compared to the theoretical nebu- 
lar continuum flux computed with the physical conditions rele- 
vant to the object. After applying the appropriate reddening the 
model flux is scaled by the ratio of observed to theoretical H^ 
line flux. Finally, the theoretical continuum is utilized for the 
nebular subtraction. Due to the good agreement between theory 
and observation, a possible stellar contribution to the observed 
Balmer jump can be excluded. 

4.2. Spectral fit 

The central star of SMP61 shows a typical [WC5]-type spec- 
trum dominated by emission lines of He, C, and O. Apart 
from a narrower line width, the spectral appearance of SMP 61 
is similar to the well-known massive WC 5 star WR 111. 
Consequently, the fit criteria and the fitting tec hnique are sim- 
ilar to the analysis of WR 111 described by ^ Grafener et alJ 
(|2002). An important difference between both objects concerns 
the low iron abundance of SMP 61. The iron forest around 
1300 A - a combination of a large number of iron emission 
and absorption lines - is barely visible in the HST observations 
of this object. 

The stellar analysis is based on fitting the absolute flux dis- 
tribution including the observed wind emission lines of He ii. 
Cm, Civ, Oiii, Oiv, and Ov. A hydrogen-free WC surface 
composition with a metallicity of 0.3 solar is adopted (refer- 
ring to the compilation of solar element abundances in iGravl 




Fig. 3. Spectral fit of the central star of SMP 61 . The observation (thin line) is shown together with the synthetic spectrum (thick 
line, grey). The model parameters are compiled in Table|3] Prominent spectral lines are identified. After correction for interstellar 
extinction, the observed and the model flux are both divided by the same model continuum, i.e. absolute fluxes are compared. 
Additionally, a correction for interstellar Ly a absorption is applied to the model spectrum. 



ill 9921 p. 3 18)). This value matches the derived nebular oxy- 
gen abundance. The iron mass fraction is correspondingly set 
to Xpe - 4 10""^, and the silicon abundance to Xsi - 2 10"^. 
The other model parameters are kept free throughout the fit- 
ting procedure. The hydrogen abundance is set to zero. This 
assumption cannot be checked empirically, because the strong 
nebular Balmer emission masks any possible stellar contribu- 
tion. However, in previous works no hydrogen is found on the 
surface of early-type [WC] stars, and even f or late types only 
small amounts of Xh < 0.1 are detected jKoesterkel l200 it 
iDe Marco & Barlowl20nil). 

The final model fit is presented in Fig.|3together with the 
derived stellar parameters in Table|3] As expected, the stel- 
lar temperature, transformed radius, and surfa ce composition 
are similar to the values obtained for WR 111 dOrafener et al.l 
while the luminosity is lower by 1 .5 dex. A good fit 



quality is obtained with some exceptions like Cm 2297 A (a 
transition that is mainly fed by dielectronic recombination), 
the classification line Cm 5696 A, or the Heii/Civ blend at 
54 12/5470 A. The latter reacts very sensitively to changes of 
the ratio of carbon to helium abundance. Additionally, our 
model with standard iron abundance produces too strong iron 
features in the UV. 

Test calculations for different iron abundances are pre- 
sented in Fig.|4] For a mass fraction around 5 10"^ (~ 0.04 so- 
lar) the overall strength of the iron forest is well reproduced, 
but the wavelengths of the single transitions still do not match 
the observation. For this reason we can only set an upper limit 
ofXpe < 1 10^4 0.07 solar) to the iron abundance. Compared 
to the standard LMC metallicity found in the nebula the central 
star shows a significant iron deficiency. In addition, the absence 
of Niv fine emission at 1721 A, which is relatively common 



6 



G. Stasinska et al.: The planetary nebula LMC-SMP 61 and its [WC]-type central star 



Table 3. Model parameters for the central star of SMP61. 
Stellar core parameters: Stellar luminosity L*, temperature 
and radius Wind parameters: mass-loss rate M, terminal 
wind velocity v^o, clumping factor D, and the corresponding 
transformed radius R^. Atmospheric abundances: Mass frac- 
tions Xhb, ^c, and Xq are determined by the model fit (Fig.|3 
whereas Xsi is set to a standard LMC value. For Xpe and Xn 
upper limits are given. The hydrogen mass fraction could 
not be determined in the present analysis. Interstellar parame- 
ters: Distance modulus M - m for the LMC, color excess due to 
interstellar extinction Eb-v, and Ly a hydrogen column density 
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Xlk 

Fig. 4. Model calculations for different iron abundances. The 
strength of the iron forest indicates a significant iron deficiency 
on the central stars surface (Xpe < 1 IQ-"*). 

among galactic [WC] stars, allows for the determination of an 
upper limit for the nitrogen surface abundance. Calculations 
with a full nitrogen model atom imply a value of Xn < 5 10"^. 

For the subsequent nebular analysis, it is important to dis- 
cuss the reliability of the derived ionizing fluxes. The main dif- 
ficulty in the present analysis arises from the nebular contin- 
uum dominating the observed flux distribution from k 2500A 
to the visible wavelength range. At first sight, the situation 
seems to be well established, because the subtraction of the 
nebular continuum works without difficulties. In addition, the 
resulting stellar flux distribution is very well reproduced by our 



models, and the extinction parameter determined from the neb- 
ular spectrum is perfectly compensating the interstellar 2200 A 
feature (see Fig.|5}- However, slight uncertainties of the flux 
level in the optical wavelength range remain. For this reason, 
the determination of the stellar luminosity primarily relies on 
the UV, i.e. it becomes extremely sensitive to the applied red- 
dening parameter /iB-y. 

Moreover, already small changes of the model parameters 
may significantly affect the derived values for the stellar lumi- 
nosity and ionizing fluxes, because the observed flux distribu- 
tion is still far away from the flux maximum at ^ 350 A. Test 
calculations indicate that a variation of the model parameters, 
leading towards an inferior but still acceptable fit quality, can 
change the derived luminosity up to 0. 1-0.2 dex. Including the 
observational uncertainties mentioned above, the possible error 
may reach values up to 0.3 dex. 

In Fig.lsj the overall flux distribution as obtained from our 
model calculations is shown together with the observed flux 
distribution, and the ionization edges of O ii. He i and H i. The 
importance of testing the ionizing fluxes for a reliable deter- 
mination of the overall flux distribution of the central star is 
obvious. A detailed analysis of the nebular spectrum is there- 
fore capable to put substantial constraints on our atmosphere 
calculations. 

5. Nebular modelling technique 

5.1. The photoionization code 

The photoionization models for the nebula were constructed 
using the code PHOTO which computes the radiative trans- 
fer of the ionizing photons in spherical symmetry with the 
"outward only" assumption for the diffuse radiation field and 
solves at each radius the equations of ionization equilibrium 
and thermal balance. The atomic data are the same as used 
by Stasiiiska & Leitherer (1996). The only modification is that 
low-temperature dielectronic recombination coefficients for S 
and Ar ions were introduced (assuming that they are equal 
to the computed values for the corresponding second-row el- 
ements, O and Ne respectively). The effects of dust are treated 
as described in Stasiiiska & Szczerba (1999, 2001). 

5.2. The Strategy 

The aim is to reproduce all the observational constraints within 
a certain degree of tolerance (to be clarified later). The observa- 
tional constraints are the emission line fluxes in the optical and 
the ultraviolet, the fluxes in the 12/vm, 25jum and 60yum IRAS 
bands and the surface brightness distribution as measured by 
HST. Out of the zone affected by the central star, the surface 
brightness distribution is dominated by the [O m] /14959, 5007, 
H/3 and Ha lines (with the [Oiii] lines giving the major contri- 
bution except, perhaps, in the outer parts of the nebula). 

The models are computed using as an input the radiation 
field from the stellar atmospheres described in the previous sec- 
tion. The other input parameters are the inner nebular radius, 
the gas density distribution, the abundances of the elements and 
the characteristics of the dust grain population. 
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Fig. 5. Overall flux distribution and 
ionizing fluxes of SMP61. The model 
flux is shown in grey together with the 
extinction-corrected observed stellar flux 
of SMP61 (black, dashed). The edge fre- 
quencies for the ionization of Hi, Hei 
and O n are indicated by vertical lines. 



We assume that dust is composed of graphite grains with a 
classical power-law size distribution of exponent -3.5 between 
0.003//m and 0.3/im. The total dust-to-gas mass ratio is cho- 
sen so as to reproduce the observed flux in the 25pm IRAS 
band (which includes nebular fine structure Unes). As a matter 
of fact, in SMP61, the observed IRAS fluxes correspond to a 
small amount of dust, which then does not aff"ect significantly 
the nebular Ylfi luminosity or its ionization structure. We do 
not introduce an additional population of small organic grains, 
as in Stasinska & Szczerba (2001), since the eff'ects of those 
grains are essentially to heat the gas and, as will be clear later, 
our initial models rather show an excess of heating. 

We first determine the density distribution that reproduces 
the observed surface brightness distribution. By trial and er- 
ror, we find that the observed surface brightness profile is well 
reproduced by a spherical model with a hydrogen density dis- 
tribution described by 



n = no exp 



1.8410" ] 



(1) 



where r is the distance to the star in cm and Hp is the density at 
the inner boundary. 

The chemical composition with which we start is the one 
derived by classical empirical methods in Sect. 2.2. (the abun- 
dances of Mg, Si, Fe with respect to H are fixed at an arbi- 
trary value of 10"*" and do not affect our results). In a first step, 
we aim at reproducing the main characteristics of the nebula, 
which can be summarized by the following quantities: F(HyS) 
the total nebular flux in HyS, the nebular outer angular radius, 
6, defined as the radius of the nebular image at zero intensity, 
and several line ratios: [Sn]/I6731/6717 to test the gas den- 
sity, [Ora] /14363/5007 to test flie gas temperature in the 0++ 
zone, [On] /l3726,3729/[On] ^7320, 7330, to test the density 
and temperature in the 0+ zone; [Oiii] .15007/[On] ^727 to 
test the ionization structure; [Oiii] /15007/H/3 to test the heat- 
ing power of the central star, and [O i] /16300/Hj8 to test the 
importance of a warm neutral zone. It is only when we are rea- 
sonably satisfied with these constraints (which we wiU refer to 



as the restricted set) that we modify the elemental abundances, 
if needed, to obtain an acceptable model. 

In order to easily visualize whether a model is acceptable, 
we compute for each constraint P the quantity 



o-p = (logPmod - log/'obs)/?, 



(2) 



where f mod is the value returned by the model, fobs is the 
observed value, and / the accepted tolerance in dex for this 
constraint. We determine the values of t for each constraint 
from the estimated uncertainties in line fluxes, including the 
uncertainty in the dereddening process. Note that we use the 
outer nebular radius as a convenient way to visualize the 
constraint on the nebular size (but in fact we require the entire 
brightness profile to be reproduced by the model). Rather 
than using line fluxes with respect to H/3, as is done usually, 
we prefer to use line ratios that are easier to interpret. Apart 
from the traditional temperature and density diagnostics men- 
tioned above, we also use [Sii] i4069,4076/[Sii] ^6717,6731 
and O III] /I1661/[0iii] /15007, which increase with nebu- 
lar temperature and density, and Cm] /I1909/Cii/14267, 
which increases with nebular temperature. The line ra- 
tios [Oiii] ^5007/ [On] ^3727, Cm] /11909/Cii] i2326, 
Nm] il750/[Nii] i6584, [S iii] ^6312/[Sii] ^6717,6731, 
and [Ar rv] ^4740/[Ariii] A1U6 
structure of the nebula. Finally, 
[Sii]/l6725/[Oii] ^3727, 
[Neiii] ^3869/ [Oiii] ^5007 and 



reflect the ionization 
[N II] /i6584/[Oii] /13727, 
Ciii]il909/[Oiii] A5Q01, 
[Arm] .l7136/[Om]^5007 



are closely linked to the abundances of the elements with 
respect to oxygen while the IRAS fluxes reflect the abundance 
and temperature of dust. The entire list of constraints used in 
this study, together with the associated tolerances is given in 
Table 4. Note that we do not consider redundant lines, such as 
[N n] A6548 whose ratio with [N ii] /16584 is independent of the 
physical conditions in the nebula. The results of our models are 
presented in graphical form in the upcoming figures, where the 
values of crp are plotted for all the constraints P. The quality of 
a model is then easily judged by the location of the points with 
respect to the line of ordinate 0. Note that we consider a model 
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as fully satisfying only if each of the values of crp is roughly 
found between -1 and +1. Any strong deviation requires an 
explanation. 



Table 4. Observational constraints and tolerances for the plan- 
etary nebula SMP 61 . 



FfHiSVlO"" 


3.71 


0.18 


e 


3.50 10"' 


0.06 


He I 5876/H/3 


1.60 10"' 


0.06 


[Sii] A613l/6in 


1.81 


0.11 


[0 III] /i4363 /5007 


7.15 10"' 


0.08 


[Nii]/l5755/6584 


2.22 10"^ 


0.30 


[O n] /17320, 7330/[O n] /13726, 3729 


1.79 10"^ 


0.08 


rSnl/14069 4076/rSnl/l6717 6731 


4.92 10"' 


0.18 


Om] /11661/[0 III] /i5007 


9.18 10"' 


0.30 


C III] /11909/C II AA261 


1.7610^ 


0.15 


[0 III] /l5007/[O n] Xilll 


9.26 


0.08 


C m] /11909/C n] /12326 


5.10 


0.11 


Nm]il750/[Nii] .16584 


2.86 10"' 


0.30 


[Siii] .i6312/[Sii] .16717, 6731 


2.03 10"' 


0.15 


[Any] A 10-'4741/[Arin]i7135 


4.55 10"2 


0.30 


[Oiii]i5007/H8 


7.41 


0.04 


[Oi]i6300/H8 


4.20 10"^ 


0.15 


[Nil] .16584/[0 II] .13727 


3.94 10"' 


0.08 


[S ii].i6725/[0 II] .13727 


7.37 10"2 


0.15 


Ciii].11909/[Oiii] .15007 


4.26 10"' 


0.11 


[Nem] i3869/[0 m] i5007 


5.67 10-2 


0.06 


[Arin]i7135/[Om]A5007 


1.1910-2 


0.10 


F(12y[im) 


8.00 10"2 


0.18 


F(25jum) 


1.3010"' 


0.10 


F(60//m) 


1.60 10-1 


0.18 



6. Results from the nebular model fitting 

6.1. First models 

We use the stellar atmosphere model described above as an 
input to the photoionization code. From the point of view 
of the nebular diagnostics available, there is no difference 
between using the model atmosphere with iron mass fraction 
Xpe = 4 10"'* or Xpe = 5 10"^. The hydrogen density distribu- 
tion of the nebular model is given by Eq. (1). As will be made 
clear later on, one of the main shortcomings of the models is 
that the excitation, as measured by [Om] /l5007/[Oii] /13727 is 
too high. We therefore discuss only models tailored to provide 
the lowest possible [Oiii] /l5007/[On] 13121 ratio. For exam- 
ple, there is no point to present models in which the nebula 
would have an inner radius smaller than 5 10*^ cm (which 
corresponds to 0.065"at the distance of the LMC). Also, while 
planetary nebulae are not necessarily ionization bounded, op- 
tically thin models are not discussed, since they lead to higher 
[Om] /15007/ [On] .13727 ratios than ionization-bounded 
models. We tried various values of the irmer density no- Note 
that the [Sii] /16731/6717 ratio is not very sensitive to no, 
due to the fact that the hydrogen density decreases outwards, 
and that, in addition, the electron density in the [Sn]/16725 
emitting zone is smaller than the hydrogen density. We show 
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Fig. 6. Representation of a series of photoionization mod- 
els of S1VIP61. The plotted quantities are the values of crp - 
(logPmod - log/^obs)/? as a function of the observational con- 
straint P (see text in Sect. 5.2). The photoionization models 
have a density structure that reproduces the observed surface 
brightness profile and they use as an input the stellar atmo- 
sphere flux described in Sect. 4. The chemical composition of 
these models is described in Sect. 2.2. All the nebular mod- 
els in this figure have a volume filhng factor of one. Models 
with inner density no=1.0 1.05 10^ 1.1 lO"*, 1.15 W cm'^ 
are represented by circles, squares, triangles and diamonds re- 
spectively. In principle, a completely satisfactory model should 
have all its points located between the dashed lines. 



in Fig. 6 a series of models with no equal to 1.010* cm"^ 
(circles), 1.05 10"* cm"^ (squares), 1.1 10"* cm"^ (triangles), 
1.15 10"* cm"^ (diamonds). AH these models have a dust-to-gas 
mass ratio of 3 10"*. We see that a value of no between about 
1.02 and 1.12 10* cm"^ leads to an angular radius (and a sur- 
face brightness profile) compatible with the observations. But 
they give too large F{WP) (by a factor about 3 compared to the 
observed value after correction for extinction). We cannot help 
by adding more grains to absorb part of the Lyman continuum 
radiation, since the modelled flux in the 25/an IRAS band is 
already larger than observed. Thus, we conclude that, either 
the real number of stellar photons in the Lyman continuum 
is smaller than in the model atmosphere or that there is some 
leakage of ionizing photons. As seen in Fig. 6, the photoion- 
ization models also predict too large [Om] /l5007/[Oii] /13727 
and [Oiii]i5007/H/3 (with a crp of 3-4), and sHghdy too 
large [Om] .14363/5007. This cannot be improved by chang- 
ing the oxygen abundance or by fine tuning the C/0 and 
S/O ratios to fit the observed Cm] .11909/ [OiiiJ i5007 
and [Sii]/l6725/[Oii] /13727 ratios. The models also 
have too low [On] .17320, 7330/[On] i3726, 3729 and 
[Sn] .14069, 4076/[S II] .16717,6731 ratios, the agreement with 
the observations being better for higher values of no, because 
the nebula is then more efficient in absorbing photons and the 
ionization front occurs at a higher value of the density. 

In order to get away with the F(HyS) problem, one 
could argue that the covering factor of the nebula is smaller 
than unity, so that some ionizing photons escape with- 
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out interacting with the gas. But this would not help in 
solving the [Oiii] ^5007 /[On] ^3727, [Oiii] ^5007/H/3 and 
[Om] /14363/5007 problems. 

We do not attempt two-sector models of different densities 
(such as those of Clegg et al. 1987 or Dudziak et al. 2000) for 
the following reasons. First, as mentioned in Sect. 2, the neb- 
ular image is extremely round. This, a priori does not mean 
that the nebula is spherical. Indeed, it could be elongated but 
seen pole on. But, in that case, the nebular extension along 
the line of sight would be larger than the one observed on 
the images. Replacing a portion of the model nebula with two 
sectors of lower density (to obtain a larger dimension) would 
considerably worsen the [On] /i7320, 7330/[Oii] /13726, 3729 
and [Sii] ^4069,4076/[Sii] ^6717,6731 ratios while not help- 
ing much for the [O iii] /l5007/[O ii] /13727 ratio and worsening 
the [O m] A5QQ1 /H/3 ratio, as we can judge from spherical mod- 
els built at lower densities. Besides, if the low density sector is 
density bounded, this may indeed help with the H/S flux, but 
will worsen the [Om] /l5007/[Oii] A3121 ratio. 

One way to lower the excitation of the nebula is to assume 
that it is clumpy. Models that differ in density and filling fac- 
tor but have the same value of the product n^e, where e is the 
volume filling factor, will have about the same external ra- 
dius but will have lower excitation for higher densities. Such 
a series of models is shown in Fig. 7, where the values of the 
couple (no, e) are (1.05 10"* cm-\ 1.0) (circles), (1.5 10"* cm^^ 
0.49) (squares), (2.0 10"* cm-^ 0.27) (triangles), (3.0 10'* cm^^ 
0.122) (diamonds). As seen in the figure, models with higher 
densities considerably improve many of the line ratios, and 
especially the [Om] /15007/[Oii] /13727 and [Om] /15007/H/3 
ratios as well as the [On] ^7320, 7330/[Oii] /13726, 3729 
and [Sn]^4069,4076/[Sii]/16717,6731 ratios and also the 
[O i] /16300/Hy6 ratio. It is then easy to fine tune the N, S and C 
abundances. A model with no between 2.0 10'* and 3.0 lO'* cm"^ 
would then be in agreement with all the observational con- 
straints (except C m] /il909/Cn/i4267, but see later) if one as- 
sumed a covering factor of roughly 1/3. 

The predicted emission in the IRAS bands deserves some 
comments. First, we note in Fig. 7 that models with lower exci- 
tation show lower fluxes in the IRAS bands. The reason is that 
the larger proportion of neutral hydrogen particles in the nebula 
reduces the amount of Lyman continuum photons absorbed by 
dust. Another point to make is that, in the hypothesis of incom- 
plete coverage of the ionizing source by the nebula (and if dust 
is intimately mixed with the gas), the predicted IRAS fluxes 
should be larger than observed by roughly the same amount 
as F(li/3) is larger than observed in our models computed for 
a covering factor of unity. Obviously, this is not the case of 
our predicted flux in the 12yum band (the observed flux in the 
60//m band is only an upper limit so there is no problem there). 
Increasing the dust abundance only worsens the problem, since 
then the dust particles are more numerous to share the photons 
that heat them, and the temperature of the dust grains becomes 
smaller, decreasing the flux in the 12yum band with respect to 
that in the 25/vm band. Actually, there may be a simple explana- 
tion to our under-prediction of the flux in the 12//m band. Our 
code does not handle PAHs, while PAHs are known to exist in 
planetary nebulae with Wolf-Rayet central stars (Szczerba et 
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Fig. 7. 



Same as Fig. 6 but for models with different filling fac- 
tors. Models defined by (no, e) equal to (1.05 10"* cm""*, 1.0) 
are represented by circles, (1.5 10^ cm"-', 0.49) by squares, 
(2.0 10'* cm-\ 0.27 by triangles, (3.0 10"* cm^^ 0.122) by di- 
amonds. 



al. 2001, Cohen 2001) and they contribute to the emission in 
the 12yum band. Note also that our predictions concerning dust 
emission are given here only for the ionized part of the neb- 
ula. The presence of dust in the neutral zone would increase 
the fluxes in the 25yL(m and 60;um bands. We do not consider 
this in the modelling presented here, since it brings no useful 
information in the problem we are interested in. In the follow- 
ing, we therefore simply require from our models that they do 
not over-predict the 12/im and 25/im fluxes with respect to the 
observations. 

We have found that with a model having an inner density 
of about 2-3 10"* cm"-' and a filling factor about 0.3 - 0.1, 
we are able to reproduce many of the observational con- 
straints, provided that we assume a covering factor of about 1/3. 
However, this model is still not entirely satisfactory. Mainly, 
the [O m] temperature is still too high and so is the predicted 
[Om] A5QQ1 /H/3 ratio. And a conspicuous shortcoming of the 
model is that the Cn/i4267/Cm] /1 1909 ratio is predicted too 
large by a factor of about 5. Note that, while the density struc- 
ture of planetary nebulae is known generally not to be uniform 
on a small scale from high resolution images, the clumps or fil- 
aments of gas are actually embedded in a more diffuse medium. 
This diffuse medium is expected to be more highly ionized, and 
thus its presence will enhance the [Om] /l5007/[Oii] /13727 
predicted by our simple model. We do not attempt to construct 
a two-phase model but we keep this in mind. 

In the following, we look for more satisfactory solutions by 
releasing some of the assumptions of our models. 



6.2. Models with modified stellar fluxes 

Due to the compact and spherical appearance of the nebula, 
the covering factor of 1/3 which is needed to explain F{li/3) 
seems to be very unrealistic. The obvious next step is there- 
fore to assume that the discrepancy is at least partly due to an 
overestimation of the stellar ionizing radiation field in Sect.|3 
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Fig. 8. Same as Fig. 7 but for models with in which the stellar 
photon flux has been divided by a factor 2 with respect to the 
original model for all the energies larger than 13.6 eV. Models 
defined by («(), e) equal to (7.5 10^ cm"^, 1.0) are represented 
by circles, (1 lO'* cm"^, 0.6) by squares, (1.2 10"^ cm"-', 0.4) by 
triangles, (1.4, 10"* cm"-', 0.3) by diamonds. 
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Fig. 9. The same as Fig. 8 for models ionized by a soft- 
ened stellar flux (see text). Models defined by («o, e) equal to 
(7.5 10^ cm"^, 1.0) are represented by circles, (1 10"^ cm"-', 0.6) 
by squares, (1.2 lO'* cm"^ 0.4 by triangles, (1.4, 10'* cm"', 0.3) 
by diamonds. 



As explained in Sect. l4.2l the eiTor for the derived stellar lumi- 
nosity may be as high as a factor of 2 under very pessimistic 
assumptions. 

Let us therefore first consider the case where we assume 
that the number of photons emitted by the star at any energy 
larger than 13.6 eV is only 50% of what is predicted by the 
model atmosphere (from the point of view of the photoioniza- 
tion model, it is equivalent to assume that the star luminosity is 
lower than given in Sect. 4 by a factor 2). The results are shown 
in Fig. 8. Here, the values of the couple no, e are 7.5 10^ cm"^, 
1.0 (circles), 1.0 10"* cm"^ 0.6 (squares), 1.2 10^ cm"^ 0.4 
(triangles), 1.410** cm"^, 0.3 (diamonds) while the dust-to- 
gas mass ratio is still assumed to be 3 lO"^"*. Since the num- 
ber of ionizing photons is smaller than in the previous case, 
we had to take smaller densities to match the observed 
size of the ionized nebula. Globally, the ionization parame- 
ter is nevertheless smaller, so that [Oiii] /l5007/[Oii] /i3727 
can now be fitted even with a model with a filling factor 
of unity. Of course, the [On] ^7320, 7330/[Oii] /13726, 3729 
and [Sii]/l4069,4076/[Sii] ^6717,6731 ratios are better re- 
produced for the highest density models of this series, i.e. no 
around 1 .5 10^ cm"-', which again require a small filling factor 
The best models of this series are at least as good as the best 
models shown in Fig. 7 without requiring a small covering fac- 
tor. Note that the [O m] /i4363 /5007 ratio is still somewhat high 
in this series of models, and adjusting the S, C and Ne abun- 
dances (by lowering them) can only make the problem worse. 

Another option to consider is a radiation field softened 
in the Lyman continuum. This would not only lower the 
[Om] /l5007/[On] /13727 ratio, but also reduce the heating. 
Fig. 9 shows a series of models identical to those of Fig. 8, but 
this time keeping 60% of the photons between 13.6 and 24.6 eV 
and 40% of the photons with energies larger than 24.6 eV. 
Clearly, the models from this series represented by the trian- 
gles are quite satisfactory also as regards the electron temper- 



ature and the [Oiii] .15007/Hy3 ratio. [S n]/l6725/ [On] .13727 
can easily be matched by reducing the original sulfur abun- 
dance by a factor 1.3 - 2). The observed fluxes in the IRAS 
25 lira, and 60 /im bands are compatible with dust-to-gas mass 
ratios of up to 3 10"^. The quality of the fit of the optical and 
UV data is not affected by the amount of dust in the allowed 
range. The only important mismatch of the model with the ob- 
servations is the C m] /11909/C ii /14267 whose computed value 
is much higher than observed, by a factor of almost 4. This 
problem is discussed in the following subsection. 

6.3. Models with inhomogeneous nebular composition 

It is a known fact that recombination lines of C, and also of 
other elements like N, O, Ne generally indicate much higher 
abundances of the parent ions than coUisionally excited lines 
(see e.g. Liu 2002 or Esteban 2002 and references therein). 
Many explanations have been proposed (see e.g. Stasiiiska 
2002 for a review). The generally prefeiTed one is to assume 
that the chemical composition of the nebulae is not uniform 
(e.g. Torres-Peimbert, Peimbert & Pena 1990, Liu et al. 2000). 
And indeed, there is direct evidence at least in two nebulae 
(A 30 and A 78) of the presence of clumps of matter enriched 
in C (Jacoby & Ford 1983, Harrington & Feibelman 1984) and 
also N, O and Ne for A30 (Wesson, Liu & Barlow 2003). The 
central star of SMP 61 being a Wolf-Rayet, this nebula is a good 
candidate for carbon enrichment. If material seen in the stellar 
wind ends up in the nebula as cool clumps, these clumps will 
strongly emit in the Cn/14267 line. We have mimicked such 
a situation by assuming that the inner zone of the nebula con- 
sists of helium- and carbon-rich material, maintaining the same 
density distribution as before. The number of combinations to 
explore is obviously very high, and we have run many simu- 
lations. In Fig. 10, we present the result of a series of simula- 
tions with no=1.4 10^ cm""*, e=0.3 (the best combination from 
Fig. 9) with the following abundances of He and C (relative to 
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Fig. 10. Models ionized by a softened stellar flux (as in Fig. 9), 
with an inner density no-lAlQ'^ cm^^and a nebular filling 
factor 6=0.3 (the best combination of Fig. 9). Now, the inner 
3 10"^ Mq of the nebula are He- and C-rich. The abundances 
of He and C (by number relative to H) are: 1.04 and 0.4 (cir- 
cles), 0.78 and 0.3 (squares), 0.52 and 0.2 (triangles), 0.26 and 
0.1 (diamonds). Otherwise, the abundance ratios with respect 
to hydrogen are (by number): H : He : C : N : O : Ne : S : Ar 
= 1: .09: 4.0 lO""* : 3.95 10"** : 2.65 lO""* : 4.26 10"^ : 4.0 10"'' : 
8 10-^. 



H by number) in the inner 3 10~"'Mq of the nebula : 1.04 and 
0.4 (circles), 0.78 and 0.3 (squares), 0.52 and 0.2 (triangles), 
0.26 and 0.1 (diamonds). These ratios correspond to the mass 
ratio of He and C in the atmosphere of the central star, as de- 
rived in Sect. 4. Otherwise, the number abundance ratios with 
respect to hydrogen are: H : He : C : N : O : Ne : S : Ar = 1 : .09 
: 4.0 10""* : 3.95 10"^ : 2.65 lO""* : 4.26 10"^ : 4.0 10"'' : 8 10"^. 
We see that the last combination provides a very good fit to the 
data, including the Cii] /I1909/Cn/14267 line ratio. Probably, 
other combinations can be found that provide a similar fit. 

7. Discussion 

7.1. Stellar ionizing fluxes 

The nebular analysis of SMP61 indicates two possible prob- 
lems concerning the theoretical flux distribution of the cen- 
tral star: The absolute number of Lyman photons is probably 
too large by a factor of 2, and the Lyman continuum radiation 
field is possibly too hard. The first problem may be discussed 
away by assuming a very pessimistic value for the uncertainty 
of the derived luminosity. However, the second point concerns 
the slope of the energy distribution in the flux maximum and 
cannot be easily resolved. 

A simultaneous solution of both problems would be 
achieved by decreasing the effective temperature Tetf of the 
model atmosphere (i.e. increasing the radius where trqss = 
2/3). In this way, photons would be distributed from the flux 
maximum to the observed wavelength range. In effect, the de- 
rived luminosity would be decreased, and a softer flux distri- 
bution would be obtained in the flux maximum. For a corre- 
sponding decrease of T^s, the stellar core temperature must be 



lowered significantly, or the mass-loss rate must be increased. 
Of course, both operations are strongly limited by the necessity 
to fit the observed spectrum of the central star. Test calculations 
show that the simultaneous fitting of the observed spectrum and 
the nebular ionizing fluxes is not possible with the model atmo- 
spheres applied in present work. 

A possible solution concerns the treatment of dielectronic 
recombination in our stellar atmosphere code. In the present 
version dielectronic transitions are treated as optically thin, and 
their contribution to the rate equations is accounted for by the 
approach of Mihalas (1971). However, when Civ recombines 
to C III in the outer part of the WC atmosphere, the ionization 
edge of C iii becomes optically thick. In this case also the cor- 
responding dielectronic transitions become optically thick, and 
their recombination efficiency is reduced because recombina- 
tion photons cannot escape. Indeed, first tests with an improved 
treatment show that the recombination from C iv to C iii is sup- 
pressed by this effect, leading to higher derived mass-loss rates 
and slightly lower stellar temperatures. A detailed investigation 
of this topic will be presented in a forthcoming paper 

7.2. The effect of He-rich clumps on the nebular 
spectrum 

He-rich clumps, which we have invoked in Sect. 6.3, can in 
fact soften the radiation field available to the rest of the neb- 
ula under certain conditions, by selectively absorbing the pho- 
tons above 24.6 eV. The total mass of helium M^^ required 
to absorb all the He°-ionizing photons emitted by the central 
star is of the order lO^'/weMo (where We is the electron den- 
sity in cm""*). If the integrated mass of helium in the clumps 
is at least equal to this value and if the clumps are located 
close to the inner boundary of the nebula with an integrated 
covering factor of unity, then the He"-ionizing photons would 
be completely blocked by the He-rich material. Such an ex- 
treme situation does not correspond to the case of SMP61, 
since [O iii] /15007 is emitted not only in the most central part 
but in the entire nebula. One can however imagine a less ex- 
treme situation where the high density He- and C-rich clumps 
are distributed over the entire volume of the nebula and soften 
the average radiation field available to the bulk of the nebular 
material. In such a case, a softening of the predicted stellar en- 
ergy distribution in the Lyman continuum would perhaps not 
be necessary. Note however that such clumps are not efficient 
in absorbing radiation between 13.6 and 24.6 eV, unless they 
are very dusty, which for SMP61 seems excluded by the ob- 
served IRAS fluxes. Obviously, the impact of these He- and C- 
rich clumps on the global spectrum of the nebula will strongly 
depend on their density, on the amount of mixing with nebular 
material and on their spatial distribution in the nebula. 

In the case of SMP 61, we can estimate the maximum mass 
of helium contained in the clumps by considering that the He/C 
proportion is the same as in the stellar wind and by using the 
observed flux in the Cii/14267 line. The observed flux in this 
line corresponds to a carbon mass Mc of about lO/weMo, if 
most of the carbon is in C^^ form and if the line is emitted at 
a temperature of 8000 K, which is the case for the C-rich zone 



12 



G. Stasinska et al.: The planetary nebula LMC-SMP 61 and its [WC]-type central star 



in our models. Since helium and carbon are roughly equal by 
mass in the stellar wind, this means that the total helium mass is 
much smaller than to M"^^, implying that there is no significant 
blocking of He*^ -ionizing photons by the clumps. If the clumps 
are of much higher density than in our models, C^^ may be 
partly recombined in the clumps, increasing Mc with respect 
to M^^, but in that case the clumps will have a smaller cross 
section and will be less efficient to block the stellar radiation. 
Therefore, in the specific case of SMP61, we do not think that 
the He-rich material will significantly soften the stellar radia- 
tion available to the nebula. 

7.3. Stellar parameters and evolutionary status 

The derived stellar parameters are in line with recent analy- 
ses of galactic [W CI-type ce ntral stars based on line -blanketed 
mode ls LDe Marco fc Crowther. .1998. 1999: De Marco et alJ 
l200lt ICrowther et alJ l2003l). In contrast to previous un- 
blanketed calculations ("Koeste rke & HamannI Il997 a*b*). the 
new models show a trend towards similar surface mass frac- 
tions for early and late [WC] subtypes, with Xq x 10% and 
Xc/Xue ~ 1. 

Due to the known distance to the LMC, SMP61 offers 
the rare chance to determine the luminosity of a [WC]-type 
central star. From the analysis of the stellar spectrum alone 
we infer a value of IQ^ '^Lq. The nebular analysis reveals that 
this value is probably too large by up to 0.3 dex - dependent 
on the assumed nebular covering factor. The luminosity of 
SMP61 is therefore in the range o f 10^-^ -10-^ '^Lp. For a sec- 
ond object, BD+303696, Li et al] (l2002l) derive a distance of 
1.2kp c from the angti l ar exp ansion of the nebula. For this dis- 
tance, [Crowthe^^U ( 120031) obtain a value of - lO^^L©. 
Both values are in the range that is expected for the major- 
ity of all post-AGB st ars, with a typical mass around 0.6 M© 
(see iBlockert 1995; .Herwigl200 lb ). Altogether, no evidence is 
found that [WC] central stars have different masses than H-rich 
objects although, admittedly, the number of cases with known 
Li, is small. 

An interesting hint on the evolutionary status of SMP61 
is given by its low iron and nitrogen abundances. For a stan- 
dard solar composition as given by 'GraV ("19921 p. 3 18) or 
fcirevesse & Sauval ( 1998) a ratio of lo g(Fe/0) = -1.33 
is expec ted. More recent investigations dAllende Prieto et all 
l2001albh imply a lower solar oxygen abundance and give 
log(Fe/0) = -1.19. The upper limit of Xpe < 1 10 that is 
derived in Sect.|3 translates to log(Fe/H) < -5.54 for a solar 
composition. In relation to the oxygen abundance of the nebula 
(log(0/H) = -3.58) we obtain log(Fe/0) < -1.96, i.e. iron 
is under-abundant by at least 0.63 dex in the atmosphere of the 
central star In addition, no nitrogen is detected in our analysis. 
The upper limit of Xn < 5 10"^' lies 0.9 dex below the abun- 
dance derived for the nebula. 

The detection of an iron deficiency is in line with recent re- 
sults from analyses of PG 1 159 stars, the proba ble descendants 
of [WC]-type central stars. iMiksa et alJ(l2002h set upper limits 
from to -1.5 dex solar for the iron abundance of 15 PG 1 159 
stars. Additionallv. .Werner et al.. (.2003.) find an iron depletion 



of at least -1.5 dex for the central star of A78, a [WC]-PG 1 159 
transition object. AlsofCrowther et al. (2003) find evidence for 
an iron abundance of 0.3-0.7 dex below solar for BD+30 3696 
and NGC40. However, for these two objects the ne bulae show 
oxyge n abundances significantly below solar (see IPefla et all 
2001), so that the iron deficiency may be attributed to a low 
initial metallicity. 

Such an iron deficiency is expected for material which 
has been exposed to s-process nucleosynthesis in the He- 
intershell of therrnally pulsating AGB or post-AGB stars (see 
iLugaro et alJ2003UHerwig et alJ2003h . This region is expected 
to consist of partially He-burned material, i.e. mainly He, C, O 
and no N. The observed surface composition of SMP 61 there- 
fore exactly resembles the abundance pattern that is expected 
for a He-intershell where the s-process has been active. 

For the formation of central stars with [WC] surface com- 
position, the H-rich layers above the He-in tershell must be 
rem oved or mixed with intershell material. iHerwi 3 (l2001al) 
and iBlockeJ (l200ll) demonstrated that the latter is possible 
by dredge-up of intershell material after an AGB-final thermal 
pulse (AFTP) or a late- or very late thermal pulse on the post- 
AGB (LTP, VLTP). In case of a LTP or VLTP very long evolu- 
tionary timescales on the post-AGB are expected, because the 
central star first evolves to the blue as H-burner, then back to 
the red during the thermal pulse, and again to the blue as He- 
burning [WC] star. 

The observed absence of nitrogen for SMP 61 strongly fa- 
vors the AFTP scenario because for a LTP or VLTP nitrogen 
is produced by the CNO-cycle in the H-burning phase. In the 
AFTP scenario the last thermal pulse occurs when the star is 
leaving the AGB, and the central star directly enters the post- 
AGB as a He-burner. However, to obtain a sufficiently high 
probability for the AFTP to occur at small enough H-envelope 
masses, a cou pling of m ass-loss to the thermal pulse cycle is 
required (Bloc keil200ll) . From the upper limit of < 5 10"^ 
we can derive an upper limit for the mass of the remaining H- 
envelope after the final thermal pulse. 

Under the assumption that nitrogen has nebular abundance 
in the H-envelope, and is completely destroyed in the He- 
intershell, it follows that the maximum mass of the H-envelope 
can be only 0.1 5 of the dredged-up mass from the intershell. 
iHerwid (12001 a' ) finds that the dredged-up mass is typically 
of the order of 6 lO^^M©. Consequently, the remaining H- 
envelope mass at the final thermal pulse must have been smaller 
than 9 10"^^^. This value is significantly smaller than the cor- 
respond ing envelope masses in the AFTP models of Herwig l 
(l2001ah . For these models the mass-loss on the AGB has 
been tuned in such a way that envelope masses of 3 IQ-^ and 
4 IO^^Mq were obtained after the thermal pulse. The present 
work implies even lower values, and therefore also lower H- 
abundances for the resulting [WC]-type central stars. 

Additional evidence against the LTP and VLTP scenarios 
is due to the very young appearance of the planetary nebula. 
Usin g the obser ved nebular expansion velocity of 29.3 kms ' 
( Vassiliadis et al. 1998) and angular radius from Table 4 we 
find an expansion time of abou t 3000 yr, whereas evolutionary 
tracks of'Blocker' (' 1995", '200 1') for a 0.625 Mq central star un- 
dergoing a VLTP imply a timescale of at least 8000 yr to reach 
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a temperature of 85 kK (dependent on the ti me whe re the zero 
poin t for the PN age is assumed). Moreover. iGornv & Tvlenda 
(1200 0) find general evidence against the LTP and VLTP sce- 
narios due to the similarity of the [WR]PN and non-[WR]PN 
populations. 

7.4. Stellar mass-loss history 

We know from the analysis of the central star atmosphere that 
the presently observed mass-loss rate is M = 10"^ '^ yr"' 
and that the carbon mass fraction is 0.52 (see Table 3). As men- 
tioned in Sect. l6.3l the total mass of carbon in the nebula re- 
quired to produce the observed flux in the C ii /14267 line, Mc, 
is about 10/«eMQ. In our model, the value of We in the zone 
where this line is mostly emitted is about 210^ cm"-' which 
implies Mq ~ 5 IO^'^Mq. Note that Mc is just sUghtly larger 
than the total mass of carbon in the C-rich zone, and that, if 
the density of the carbon-rich clumps is larger, the total carbon 
mass will be smaller. We thus infer that the star has spent about 
1300yr in a similar state of mass-loss. 

This timescale is significantly lower than the time expected 
for a PN central star to reach T^ff = 85 kK. Evolutionary 
track s for He-burning stars of t he corresponding luminosity 
from Vassiliadis & Wood (1994) give values around 3000 yr 
IVassiliadis et al, C1998 ) show that this timesca le fits the n ebu- 
lar age of SMP61. H-burning tracks from Rlo ck^dl 99.51) im- 
ply an age of 3000 yr for a post-AGB mass of 0.625 M© (with 
L* = 10^ ''Lo) and 5000yr for 0.605 Mo (with = lO^^Lg). 

Our modelling then suggests that the mass-loss of SMP61 
may be intermittent, which would indicate that it is initiated by 
other processes than radiative acceleration. Interestingly, our 
atmosphere models would confirm this hypothesis, because the 
force due to radiation pressure is much too low to explain the 
observed mass loss: For the massive WC star WR 111, which 
has a very similar spectral appearance to SMP 61, the radiation 
pressure, as calculated in our models, supplies abo ut one half 
of the energy necessary to drive the stellar wind (.Grafener et al.l 
I2OO2I) . On the other hand, for SMP 61 the same models provide 
only 17% of the wind energy. However, one must be aware that 
the timescales for mass-loss and evolution derived above both 
suffer from large uncertainties. The spectroscopically derived 
mass-loss rate depends on the clumping factor (M oc 1 / ^/D 
for constant Rt, see Sect. l3.2> which is only roughly known. 
The error margin for M may therefore be as large as +50%. 
The stellar evolutionary timescale is even more uncertain, be- 
cause it depends on the stellar luminosity and the mass of the 
stellar envelope above the burning shell. The latter is strongly 
dependent on model assumptions. The nebular age is also un- 
certain, mostly because the nebular velocity changes during 
the course of evolution as shown by dynamical simulations 
(Mellema 1994, Villaver et al. 2002). 

8. Summary 

We have obtained HST STIS imaging and spectroscopic obser- 
vations of the Magellanic Cloud planetary nebula SMP 6 1 and 
its central star. We have also searched the literature for other 
observational data concerning this object, in order to obtain a 



full set of observational constrains and to assess their uncer- 
tainties. 

We have performed a detailed spectral analysis of the cen- 
tral star For that, we used the Potsdam code for expanding 
atmospheres in non-LTE, which includes the statistical equa- 
tions for many ions, including ions from the iron group. The 
observed stellar spectrum is well represented by a model with 
the following parameters: = 10^ ''^Lq, = 0.42 Rq, = 
87.5 kK, M = 10-^1- Mo yr"', Voo = 1400 km s"', clumping fac- 
tor D = 4 and Ri - 4.63 Rq. The abundances of the elements by 
mass are Xee = 0.45, Xc = 0.52, < 5 Xq = 0.03, 
Xpe < 1 10 The low iron and nitrogen abundances indicate 
that the star has undergone an AGB final thermal pulse with a 
very low H-envelope mass (< 9 IO'^^Mq). The iron depletion 
is probably due to s-process nucleosynthesis during previous 
thermal pulses. 

The fluxes from the model stellar atmosphere were used 
as an input to our photoionization code to construct photoion- 
ization models of the nebula. We considered all the available 
observational constraints, within their error bars. These con- 
straints are numerous and allow a detailed discussion of the 
models. We find that the fluxes from the stellar atmospheres 
allow to reproduce many features of the nebular emission. 
However, the observed nebular properties are better explained 
if the central star produces less photons in the Lyman contin- 
uum than the model atmosphere by about a factor 2. An even 
better fit is obtained if the Lyman continuum radiation is also 
softer than predicted by the stellar model. A possible solution 
involves the treatment of dielectronic recombination in the stel- 
lar atmosphere code and wiU be presented in a forthcoming pa- 
per 

We also find that the nebula must contain extremely 
carbon rich clumps in order to reproduce the observed 
Cm] /1 1 909 /C 11/14267 ratio. The observational constraints do 
not allow us to describe these clumps in detail (density, chemi- 
cal composition, location) but these clumps are likely produced 
by ejecta from the clumpy stellar wind, possibly partly mixed 
with the nebular material. Obviously, high signal-to noise and 
high spatial resolution observations of nearby planetary nebu- 
lae excited by a similar star would be extremely useful to better 
understand the process of enrichment of planetary nebulae by 
freshly made carbon (and helium). 

From a comparison of the estimated total amount of carbon 
in the nebula with the carbon mass-loss rate measured in the 
central star, we infer that the strong WC mass-loss may have 
been active for only a limited period in the post-AGB evolu- 
tion of SMP 61. However, this last result relies on uncertain 
timescales and is subject to revision. 
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Table 2. Nebular dereddened fluxes from different authors 



lambda ion Dereddened fluxes relative to 11/3 

STIS (this work) PRT-P97a MD91 Vass98a MBC88 



1240 


Nv 


stell 










1336 


Cii 








26.0 




1406 


Si IV], O IV] 


stell 






12.9 




1550 


Civ 


stell 


stell 




72.1+50.4 




1640 


He II 


stell 






31.2 




1649 


Civ 


stell 






10.3 




1666 


[Oiii] 


6.8: 






7.4 




1750 


Nm] 


0.09 










1909 


Cm] 


316 


401 




392.7 




2297 


Cm 


stell 






30.7 




2321,26 


Cii 


62 






87.9 




2405 


Civ 


stell 






8.1 




2470 


[On] 


6.2: 






9.0 




2524 


Civ 


stell 






11.9 




3203 


He i,He n 


stell 






3.4 




3726+29 


[On] 


60 


80.0 


85.8 




<80 


3835 


H9 


7.2 


7.4 


10.4 






3869 


[Nem] 


43 


41.6 


50.6 




47.2 


3889 


Hei+H8 


19 


18.7 


24.6 






3967+70 


[Neiii]+H7 


28 


29.0 


37.9 






4026 


He i+He n 


2.6 


2.5 


3.9 






4069 


[Sn] 


2.2 


2.9 


3.9 






4102 


m 


26 


25.2 


31.2 




13.5 


4269 


Cn 


1.8 


1.4 


1.7 






4340 


Hy 


46 


46.6 


46.6 




39.9 


4363 


[Oiii] 


5.3 


5.7 


6.4 




5.5 


4471 


He I 


5.2 
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4686 


Hen 


stell 


stell 








4711+13 


[Ariv]+Hei 
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1.3 
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0.2: 
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[At IV] 




0.2: 
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100.0 


100.0 




100.0 
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He I 


1.5 


1.4 
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4959 


[Oiii] 
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222.5 




278.6 


5007 


[Oiii] 


sat? 


sat 


sat 




800.0 
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[Ni] 




0.3: 
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0.4 
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[CI III] 




0.2: 
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[CI III] 
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He I 
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18.8 




15.2 
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4.3 


6.1 




3.3 
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[Sm] 
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1.3 
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[Oi] 




1.3 


2.3 
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[Nn] 




15.8 


15.0 




7.7 


6563 
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285.8 


282.1 




282.4 
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31.5 
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22.3 
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6717 


[Sn] 
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2.2 
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3.8 
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5.3 
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He I 




8.2 


13.7 
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[Am] 




8.8 


16.3 






7320+30 


[On] 




14.3 








log F(H8) 




-12.80 


-12.61 








c(m 




0.12 


0.18 


0.19 




1.07' 


slitC) 




0.2x52 


2;10 






3.6x3.6 



' value affected by atmospheric dispersion 



